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(54) Method of analysis of samples by determination of the distribution of specif ic brightnesses 
of particles 

(57) A method for analyzing samples by measuring 
numbers of photon counts per defined time interval in a 
repetitive mode from light emitted, scattered and/or 
reflected by particles in said sample, and determining 
the distribution of the number of photon counts per said 
time intervals, characterized in that the distribution of 
specific brightnesses of said particles is determined 
from said distribution of the number of photon counts. ~ 
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Description 



The present invention relates to a method for ana- 
lyzing samples by measuring numbers of photon counts 
per defined time interval in a repetitive mode from light s. 
emitted, scattered and/or reflected by particles in the 
sample and determining the distribution of the number 
of photon counts. 

The first successful studies on fluorescence inten- 
sity fluctuations were performed by Magde, Elson and w 
Webb (Biopolymers, Vol. 13, 29-61, 1974) who demon- 
strated the possibility to detect number fluctuations of 
fluorescent molecules and established a research field 
called fluorescence correlation spectroscopy (FCS) 
FCS was primarily developed as a method for determin- rs 
mg chemical kinetic constants and diffusion coefficients. 
The experiment consists essentially in measuring the 
variation with time of the number of molecules of spe- 
cific reactants in a defined open volume of solution. The 
concentration of a reactant is measured by its fluores- 20 
cence from a small measurement volume. The meas- 
urement volume is defined by a focused laser beam 
which excites the fluorescence, and a pinhole in the 
image plane of the microscope collecting fluorescence 
Intensity of fluorescence emission fluctuates in propor- 25 
tion with the changes in the number of fluorescent mol- 
ecules as they diffuse into and out of the measurement 
volume and as they are created or eliminated by the 
chemical reactions. Technically, the direct outcome of 
an FCS experiment is the calculated autocorrelation '. 30 
function of the measured fluorescence intensity. 

An important application of FCS is determination of 
concentrations of fluorescent species having differ^ . 
diffusion rates, in their mixture. In order to separate the 
two terms in the autocorrelation function of fluores- 35 
cence intensity corresponding to translation diffusion of 
two kinds of particles, at least about a two-fold differ- 
ence in diffusion time is needed, which corresponds to 
an eight-fold mass difference of particles. Furthermore 
even if one succeeds in separating the two terms in the 40 
autocorrelation function of fluorescence intensity, it is 
yet not sufficient for determining the corresponding con- 
centrations except if one knows the relative brightness 
of the two different types of particles. 

Whereas conventional FCS yields rather limited as 
information about aggregate sizes from a simple auto- 
correlation function of fluorescence intensity, fluctua- 
tions, possible biophysical applications demand the 
ability to analyze complex mixtures of different species 
For that purpose, Palmer and Thompson studied higher so 

order correlation functions of fluorescence intensity fluc- 
tuations and have outlined methods for determining the 
number densities and relative molecular brightness of 
fluorescence of different fluorescent species (Biophys 
J- . Vol. 52,. 257-270, August 1987). Their technique may ss 

in principle proof useful in detecting and characterizing 
aggregates of fluorescent-labeled biological molecules 
such as cell surface receptors, but has a major disad- 



vantage of being rather. complex, so that data process- 
ing of an experiment including the calculation of high- 
order correlation functions lasts hours. 

A considerably less complicated method than cal- 
culation of high order auto-correlation functions for 
characterizing mixtures of fluorescent species of differ- 
ent specific brightness is calculation of higher order 
moments of fluorescence intensity out of the experi- 
mentally determined distribution of the number of pho- 
ton counts. This method was presented by Qian and 
Elson (Biophys. J., Vol. 57, 375-380, February 1990- 
Proc Natl. Acad. Sci. USA, Vol. 87, 5479-5483, July 
1990). In their demonstration experiments signal acqui- 
sition times of about 7 minutes were used for a relatively 
favourable experimental system of two kinds of fluores- 
cent particles which differed by 30 -fold in their specific 
brightness, the mixture of monomers and 30-mers. The 
method of moments is relatively simple and fastin cal- 
culations, but it allows to determine only a limited 
number of unknown parameters characterizing the 
sample because usually only about three or four first 
moments of fluorescence intensity can be calculated 
from the experiment with precision sufficient for further 
analysis. Because of this reason, the method of 
moments is hardly suitable for characterizing complex 
samples or selecting between competing' models of the 
sample or checking whether the given model is appro- 
priate. 

One object of the invention is to obtain reliable infor- 
mation /about a' sample containing particles emitting 
scattering and/or reflecting photons, which renders pos-' 
-.sible an .analysis of the sample with respect to certain 
ingredients or with respect to certain states of the sam- 
ple. 

Another object of the present invention is to sub- 
stantially extend the useful information obtainable from 
the experimentally determined distribution of the 
number of photon counts. 

The objects of the present invention are solved with 
the method having the features of claim 1 . 

What is meant by the term "specific brightness" of 
particles in the sense of the present invention is a phys- 
ical characteristic expressing in what extent a given type 
of particles is able to emit, scatter and/or reflect light. It 
is thought to characterize single particles and therefore 
the value of specific brightness is not depending on con- 
centration of the particles, neither on the presence of 
other particles. Thus, a change of the total count rate of 
photons emitted, scattered and/or reflected from the 
measurement volume, if only due to a change in con- 
centration of the total number of particles, does not 
influence the measured value of specific brightness- and 
the value of the ratio of numbers of particles of different 
types determined by the present invention. Specific 
brightness of a type of particles is usually expressed in 
terms of the mean count rate per particle which is a 
weighted average of the count rate over coordinates of 
the particle in the measurement volume. 
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The importance of the present invention for the 
analysis of samples may be illustrated by the following, 
non-limiting example: Assuming that a solution contains 
a quantity (a) of one type of particles (A) with a respec- 
tive specific brightness (la) and a quantity (b) of another 
type of particles (B) with a respective specific bright- 
ness (lb), the overall count rate of photons emitted by 
the solution depends on the expression !a*a + lb*b. 
Thus, by mere determination of the overall count rate, it 
is not possible to dissolve the value of a and/or b. Gen- 
erally, in fluorimetric measurements, the overall count 
rate of at least one type of particles is determined in an 
independent experiment. If the total number, a+b of par- 
ticles does not change with respect to this measure- 
ment, the ratio a/b or its inverse can be determined by 
mere determination of the overall count rate of the mix- 
ture in a second measurement. However, the assump- 
tion, that the total number a+b does not change 
between the two measurements, is often wrong. For 
example, adsorption effects of particles to surfaces may 
occur. Fluorimetric. measurements cannot verify the 
total number of particles a+b independently. The 
present invention overcomes these restrictions. From 
one measurement, the numbers of particles a and b can 
be determined without any prior information of their 
respective specific brightnesses. 

It is to be understood that the following description 
is intended to be illustrative and not restrictive. Many 
embodiments will be apparent to those of skill in the art 
upon reviewing the following description. By way of 
example, the invention will be described primarily with 
reference to measuring numbers of photon counts from 
light emitted by f iubrescently labelled particles in sam- 
ple. For example, in some embodiments it may be desir- 
able to measure numbers of photon counts of other 
origin than fluorescence. 

the present invention provides a method for calcu- 
lating the expected distribution of the number of photon 
counts corresponding to given real equipment and sam- 
ples of given composition. The ability to predict the dis- 
tribution of the, number of counts corresponding to 
samples of given composition allows, when studying 
samples of unknown composition, to find out the model 
of the sample yielding the closest fit between the calcu- 
lated and the experimentally determined distribution of 
the number of. photon counts. What is meant by the 
composition of the sample here is the specific bright- 
ness and concentration of fluorescent particles present 
in the sample. For example, a solution of a single fluo- 
rescent dye is characterized by two parameters: the 
concentration and specific brightness of the dye mole- 
cules. A mixture of two fluorescent dyes is characterized 
by four parameters: the concentrations and specific 
brightnesses of the two kind of molecules. A complex 
mixture can be characterized by the distribution function 
of concentration versus specific brightness of particles. 
Conveniently, concentration is expressed as the mean 
number of particles per measurement volume, and spe- 



cific brightness of particles is expressed as the mean 
count rate per particle. 

From the other side, the distribution of the number 
of photon counts depends not only on composition of 
5 the sample, but also on equipment: first of all, on the 
*- spatial brightness function characteristic to the optical 
set-up, and on some characteristics of the detectorlike 
its dark count rate, its dead time and probability of after- 
pulsing. In the interest of a high quality of analysis, 
jo which is indicated by achieving a close fit-between the 
. experimentally determined and the calculated curves, it 
is preferred to characterize the equipment adequately. 

Claim 1 does not cover the method described by 
Qian and Elson, as their method compares estimated 
75 and calculated moments of the distribution of light inten- 
sity, but' not directly the distribution of the number of 
photon counts. Qian and Elson's teaching misses the 
teaching of the present invention. 

According to the invention, a new quality of analysis 
20 of samples containing particles which emit, scatter 
and/or reflect light becomes possible. In a first step, 
numbers of photon counts from light emitted by the par- 
ticles in the sample are measured per defined time 
interval in a repetitive mode. A series of different time 
25 intervals can also be used for a more complex analysis. 
In a second step of the method according to the inven- 
tion, the distribution of the number of photon counts is 
determined, which means that it is determined how ; 
many times one has obtained a certain number of phq- 
. 30 ton counts. In a third step, the experimentally deter-; 
■ mined distribution of photon counts is analyzed directly 
' without intermediate steps to obtain the distribution of 
' 'specific brightnesses of the particles in the sample. ■ . 
In the sense of the present invention, particles are 
35 preferably luminescentiy labelled or unlabelled 'mole: 
cules or macromoiecules, or dye molecules, molecular 
aggregates, complexes, vesicles, cells, vi rusesr bacte- 
ria, beads, or mixtures thereof. 

In a further preferred embodiment, the particles 
'40 each carry : a number of binding sites for luminescent 
particles: Luminescent particles can directly or via sec- 
ondary molecules bind to these binding sites. Since 
highly luminescent particles are generated when many 
• luminescent particles bind to the binding sites of the first 
45 particles, the method according to the invention is able 
to distinguish easily between particles with a large dif- 
ference in luminescence intensity, so that even a small 
amount of bound luminescent particles can be meas- 
. ured in 'presence of an excess concentration of free 
so luminescent particles. This embodiment provides a new 
analysis of particles which do not carry a luminescent 
label by binding to a second particle which is lumines- 
centiy labelled, but whose brightness does not change 
upon binding. A commercially very important applica- 
55 tion of this method is the measurement of ftuorescently 
labelled antibodies binding to an antigen, while the anti- 
gen is binding to at least one of the multiple binding 
sites of the particle which is preferably a bead, or vice 
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versa. The method can also be applied to other types of 
interactions such as nucleic acid hybridization or pro- 
tein/nucleic acid interaction. The invention can also be 
applied for the analysis of distribution characteristics of 
said particles, such as for quality control and process 5 
control of polymers or oligomeric suspensions of parti- 
cles/ In addition, surface areas, of particles can be ana-' 
lyzed as well as distributions of surface areas of 
particles.- 

In a preferred embodiment, one. type of particle - w 
subsequently denoted A, carries more than one binding 
site. Another, luminescent type of particles, subse- 
quently denoted C, can bind (i) either directly to at least " 
one of the binding sites of particle A, or (ii) binds to at 

least one binding site of a molecule B, which in turn bind is 
to at least one of the binding sites of particle A These 
bindings may be mediated either 6y naturally occurring 
binding sites on the particles, or mediated by introduc- 
tion of specif ic binding sites to the. particles A, B and/or 
C: Since in both cases more than one of the particles of '20 
type C may bind to particle A. the complex will emit 
more photons than free particles of type C. This embod- 
iment provides a convenient way to measure binding of 
particles of type B to a particle of type C or A, although 
the particle of type B is not luminescently labelled. 25 

In a further preferred embodiment, the measure- 
ment volume is only a part of the total volume of the 
sample and said particles are diffusing and/or being 
actively transported into and out of said measurement 
volume and/or the sample is actively transported and/or 30 
optically scanned. If the fluorescent particles are suffi- 
ciently small, then diffusion is fast enough for data 

acquisition from a great number of independent meas- 
urement volumes, and data acquisition using time aver- 
aging is nearly identical to ensemble averaging 35 
However, if the characteristic time of diffusion is sub- 
stantially longer than the time interval necessary for 
measuring fluorescence intensity (which is usually 1 0 to 
50 jis), then active transport (flow or scanning) can con- 
siderably save time of data acquisition. 40 

A sample is usually characterized by values of con- 
centration and specific brightness of one or more types 
of fluorescent particles. In cases when one or more of 
these values are known beforehand, the goal of analy- 
sis is to determine unknown values, either those of con- 45 
centration, or specific brightness, or both. 

Two alternative methods for selecting the model 
yielding a fit between the experimentally, determined 
and calculated distributions of the number of photon 
counts can be used. In one embodiment the well- so- 
known least squares fitting method, where the sample is 
described by a finite (usually small) number of parame- 
ters, can be employed. The purpose is to find values of 
the parameters yielding the closest fit between the 
experimental and the calculated curves. According to 55 
the invention, values of concentrations and/or specific 
brightnesses of a number of types of fluorescent parti- 
cles can be estimated. In a further embodiment, another ' 



general method called inverse transformation with linear 
regularization (ITR) can be employed. ITR describes 
the sample using a semi-continuous distribution func- 
tion of particles versus their specific brightness and 
searches for the closest fit demanding that the. solution 
is a smooth function. (For the method of ITR, see eg 
W. H. Press, S. A. Teukolsky, W. T Vetterling. B. PFIah- 
nery Numerical recipes in C: the art of scientific com-' 
putmg. second edition,. Cambridge University Press 
1992, p. 808) ,. . . \ _ . . V 

In the following/the invention is further illustrated in 
a non-limiting manner. Particularly, it is described how 
the expected distribution of/the number of photon 
counts is calculated. , 

A first step is the calculation of the probability distri- 
bution of the number of photon counts from a spatial 
section of the measurement volume with a constant 
value of spatial brightness. This can be done by using a 
single type of fluorescent particles. It is well known that 
the probability distribution of the number of particles in 
an open volume is Poissonian, Also, if the number of 
particles inside the spatial section) is given, the number 
of detected photons per sampling time interval is Pois- 
son distributed; Consequently the overall distribution of 
the number of photon counts corresponding to the spa- 
tial section of constant brightness and an ideal detector 
is compound Poissonian. 

As the next step, one may study the case in which 
the measurement volume is divided , into a number of 
spatial sections of constant brightness. If the values of 
volumes and spatial brightnesses in each: of the sec- 
tions are known, the distribution of the number of pho- 
ton counts corresponding., to each section can be 
calculated separately All these^distributions are com- 
pound Poissonian. Furthermore, if distributions of the 
number of photon counts for all sections were known 
the overall distribution can be calculated using the fact 
that the total number of counts is the sum of the number ; 
of counts originating from different sections of the " 
• measurement volume. 

As the following step, one may study the mixture of 
types of fluorescent particles having different values of * 
specific brightness. Each spatial section of the meas- 
urement volume can be divided into a number of 
abstract subsections each containing only particles of a 
single type. A similar procedure can be applied how as 
described above for spatial sections of the measure- 
ment volume in order to calculate the overall distribution 
of the number of photon counts. 

t An experimentally, determined distribution of the 
number of photon counts is ruled not only by properties 
of the light beam, but is influenced also by nonideal 
properties of the photon detector. Statistically, the dark 
.counts of the detector behave in the same manner as 
photon counts from background light of constant inten- . 
sity. Their contribution are photon counts of Poisson dis- 
tribution. Also, the way how the dead time of the 
detector and^its afterpulsing distort the distribution of 
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photon counts are known from literature on photon sta- 
tistics (see e.g. B. Saleh, Photoelectron Statistics, 
Springer, Berlin, 1978): 

In summary, the expected distribution of the 
number of photon counts is determined, from one side, s 
by characteristics of the sample (concentrations and 
specific brightnesses of fluorescent particles of different 
kind), and, from the other side, by characteristics of the 
equipment (the sampling time interval, the spatial 
brightness function, the background count rate, the w 
dead time and the afterpulsing probability of the detec- 
tor). V 

In one embodiment, both the dead time and the 
afterpulsing probability of the detector are determined - 
from experiments in which the distribution of the number is 
of photon counts corresponding to light of constant 
- intensity is determined. The experimental parameters of 
the detector- may also be determined by measuring 
numbers of photon counts per defined time intervals in 
a repetitive mode while the detector is exposed to high 20 
frequency laser pulses. Correction for the dead time of 
the detector may be performed on the basis of a formula 
derived by Cantor and Teich (J. Opt. Soc. Am. 65, 
786,1975; see also B. Saleh, p. 272-277). Mathematics 
of afterpulsing is simple: each photon pulse can be fol- 25 
lowed by another (artificial) pulse; this happens usually 
with a constant probability. 

According to the invention, it is preferred that the 
spatial brightness function is characterized using exper- 
iments on a single type of particles. For example, if the 30 
laser wavelength 514.5 nm is used, then a solution of 
Rhodamine 6G is a convenient sample which can be 
used for characterizing the spatial brightness function! ' 

What characteristics of the spatial brightness func- 
tion can be employed when calculating the expected 35 
distribution of the number of counts are values of vol- 
umes of the sections of the measurement volume corre- 
sponding to a selected set of values of the spatial 
brightness. Typically, a Set of twenty values of the spa- 
tial brightness positioned at a constant distance from 40 
each other in the logarithmic scale have been selected, 
covering two orders of magnitude. Contribution from the' 
lower brightness areas can be accounted for by a single 
parameter, their relative contribution to fluorescence 
intensity. Intensity fluctuations of this light can be 45 
neglected. Because of the large number of the sections 
of the measurement volume, it would be less preferred 
to consider volumes corresponding to each of the sec- 
tions as independent variables. It is convenient to con- 
sider them as variables depending on a few other so 
parameters, and determine the values of these parame- 
ters which yield the closest fit between the experimen- 
tally determined and the calculated distribution of the 
number of photon counts. Conveniently, a relatively sim- 
ple model of the optical set-up is applied, which is not ss 
accounting for aberrations of the optics used, and which 
determines volumes of the sections of the measure- 
ment volume. For instance, the volumes of the sections 



depend on values of the convergence angle of the laser 
beam and the size of the pinhole. 

In fluorescence studies, it may be advantageous to 
take measures for reducing the background count rate, 
arising from Raman scattering in the solute material and 
dark count rate of the detector, with respect to the count 
rate per particle. In particular, it is in some cases pre- 
ferred to use measurement volumes smaller than 100 
pm 3 , more preferably about 1 Mm 3 . Advantageously, the 
high signal to background count rate and the small opti- 
cal measurement volume may be achieved by using at 
least one high-aperture microscope objective in a con- 
focal manner for both focussing the incident laser beam 
arid collecting light emitted, scattered and/or reflected 
by particles in said sample. 

In a preferred embodiment of the method, multiple 
photon excitation is used to excite a particle. Multiple 
photon excitation means that the sum, difference or any 
combination of wave frequences of two, three or more 
photons is used for excitation of e.g. luminescence or 
second order Raman scattering. Such an excitation 
scheme has an advantage in the sense that the excita- 
tion probability is not linearly dependent on excitation 
intensity, but on the second or higher power. Thus, the 
multiple photon excitation is mostly limited to the volume! 
of the laser focus, whereas outside the laser focus no 
spurious excitation is generated. The present invention 
profits from such an excitation scheme in the sense that 
less background is generated compared to single pho- 
ton excitation, and that there is no pinhole necessary to 
restrict the measurement volume. Thus, -the pinhole 
diameter and its imaging on the detector do not enter as 
modelling parameters. in the spatial brightness function . 
■ any more. 

In a further preferred embodiment, near field optical 
microscopy is used for focussing the excitation light of 
the particles, and/or collecting the light emitted by the 
particles. Near field optical microscopy means here that 
the light passes through an aperture with at least one of 
its dimensions being smaller than the wavelength of the 
light used and which is in direct contact to the measure- 
ment volume. The aperture may consist of an opaque - 
layer with at least one hole of said diameter or at least . 
one slit of appropriate width and/or a tapered glass fiber 
or wave guide with a tip diameter of said width, option- : 
ally coated with an opaque layer outside. 

Another preferred embodiment combines near field 
optical microscopy for the excitation light path, and con- 
ventional optical microscopy for the emission light path, 
or vice versa. The present invention profits from such a 
realization in the sense that the size of the measure- 
ment volume is reduced compared to conventional con- 
focal microscopy. Thus, the present invention can be 
used to measure higher particle concentrations as with 
other optical schemes. 

According to the invention, it may in some cases be 
preferred to select the width of the sampling time inter- 
val in such a way that in average more than one, prefer- 
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ably one to ten, photon counts per particle are yielded. < 
It may further be preferred that the width of the time 
intervaljs in average smaller than the characteristic cor- 
relation time of light intensity f luctuations. 

If more than two unknown parameters of the sam- 
ple have to be estimated, it is preferred to have no.more. 
than a few, preferably about one particle per measure- 
ment volume. 

in one embodiment, at least one individual particle 
.. is statistically analyzed in terms of its specific bright- 
ness which may fluctuate or otherwise change. 

The method of the present invention is particularly 
advantageous since not only that information losses 
and distortions are kept minimal. One new quality 
attainable by the present invention is that the data 
processing depends less from a definite mathematical 
model of the sample than do the other techniques which 
are state of the art. This means that the method is more 
reliable in terms of long term stability of an instrumental 
- realization, and that any disturbation of the measure- 20 
ment volume by e.g. turbid samples or particles inside 
the laser beam does not significantly influence experi- 
mental results. 

A further new quality attainable by the present 
invention is that the signal-to-noise ratio is much better 2s. 
compared to the techniques of the prior art. This means 
that experiments can be made within significantly 
shorter time (up to 100 fold shorter) than previously, 
showing the same signal-to-hoise ratio as previous long 
term experiments. Since photo-bleaching. of fluorescent 30 
molecules or fluorescently labelled molecules is an 
unsolved problem so far with any applied measurement 
technique in this field, especially when applied to cells, 
one is restricted to short measurement times. Thus, 
compared to the prior art, the present invention is 35 
advantageous for measurements inside cellular sys- 
tems. 

In one preferred embodiment/the present invention 
is realized in the field of fluorescence intensity fluctua- 
tion studies. The optical equipment is a conventional ,ao 
confocal FCS microscope equipped with a cw laser of 
visible light. The excitation laser beam is focussed into a 
sample which is a homogeneous water solution of a low 
concentration, typically in the nanomolar range, of fluo- 
rescent material. Fluorescence emission from a micro- 45 
scopic confocal volume of about 1 ^m 3 is collected on a 
photon detector. The measurement time which is typi- 
cally 1 to 60 seconds is divided into hundreds of thou- 
sands time intervals of typical width of 10 to 50 ps. The 
highest number of photon counts typically obtained in so 
this experimental realization of the invention herein 
described is between 1 0 and 100. 

The nature and advantages of the invention may be 
■ better understood oh the basis of the following example 
where a mixture of rhodamine dyes is analyzed. Figures 55 
1 to 7 illustrate consecutive steps of the analysis and 
their results. 

Fig. 1 . Distributions of the number of photon counts 



experimentally "determined at constant light intensities, 
time interval of 10 us and data collection time of 50 s' 
From curve fitting, the dead time of the detector was 
estimated to be 28 ± 4 ,ns; afterpulsing probability 
5 0.0032 ± 0.0008. In the lower graph, weighted residuals 
of the curve fitting are presented. • 

Fig. 2. Distribution of the number of photon counts 
experimentally determined for a solution of rhodamine 
6G at time interval of 40 us and data collection time of 
' 10 , 50 s. ' " • ' • • . • • . • • . 

Fig. 3. Normalized sizes of volumes of the twenty 
spatial sections of the measurement volume of the high- 
est brightness. Section 0 corresponds to the maximal 
value of the spatial brightness while section 19 corre- 
15 sponds to about 100 times lower brightness. Sizes of 
volumes are determined from experiments on single flu- 
orescent species. 

Fig. 4. Residuals of curve fitting corresponding to 
an experiment on rhodamine 6G solution (the experi- 
ment graphed by Fig. 2). 

Fig. 5. Distribution of the number of photon counts 
experimentally determined for three samples at time 
interval of 40 us and data collection time, of 50 s. The 
distribution corresponding to rhodamine 6G is the same 
as in Fig. 2. " . . , 

Fig. 6. The results of the inverse transformation with 
linear regularization applied to the data of Fig. 5. 

Fig. 7. Residuals corresponding to analysis of an 
experiment on the mixture, solution- of rhodamine 6G 
and tetramethylrhodamine (measured data in Fig. 5). 
Graph a: expected curve was obtained by inverse trans-' 
formation with linear regularization.< Graph b: expected 
curve was, obtained by the least squares fitting of the 
. experimental data. Graph c: residuals of the least 
squares curve fitting under a wrong assumption of sin- 
gle species. 

As the first preparatory step of analysis, the dead 
time and the afterpulsing probability of the photon 
detector are estimated. This was done by determining - 
the distribution of the number of photon counts under 
illumination of the detector by light of constant intensity. 
Since the dead time distortions are most noticeable at 
high count rates while the afterpulsing distortions are 
better resolved at low count rates, the values of the 
dead time and the probability of afterpulsingwere deter- 
mined by jointly fitting distributions of the number of 
photon counts determined at a relatively high and at a 
relatively low illumination intensity/The experimentally 
determined count number distributions are presented in 
Fig. 1, together with residuals of the curve fitting. Values 
of the estimated parameters for the photon detector 
which have been used are: dead time 28 ns, afterpuls- 
ing probability 0.003. 

The background count rate of the equipment is 
determined by measuring the count rate when the sam- 
ple holder is filled with pure water. 

As the second preparatory step, the spatial bright- 
ness distribution corresponding to a given optical set-up * 
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was characterized. For that purpose, the distribution of 
the number of photon counts corresponding to a solu- 
tion of rhodamine 6G was experimentally determined 
(Fig. 2). If the spatial brightness distribution is appropri- 
ately characterized, then the calculated curve fits the 
experimental curve. In order to numerically calculate the 
expected distribution of the number of photon counts, 
values of twenty one parameters characterizing the 
spatial profile are needed in our program: twenty sizes 
of volumes corresponding to twenty spatial sections of 
the highest values of spatial brightness, and the relative 
contribution to fluorescence light originating from areas 
of lower spatial brightness. In order to calculate values 
of these unknown parameters, a simple model of the 
optical equipment not accounting for aberrations was 
taken into use. As illustrated by Fig. 3, the determined 
sizes of the twenty volumes are reproducible, even if 
other species than rhodamine 6G are used. 

Having determined values of the twenty one param- 
eters characterizing the spatial brightness distribution in 
the way just described above, the calculated distribution 
of the number of photon counts indeed fits the experi- 
mentally curve, see Fig. 4. 

After the preparatory steps described above the 
equipment is ready for analysis. In Fig. 5, distributions of 
the number of photon counts corresponding to three dif- 
ferent samples are presented. In Fig. 6, the results of 
the inverse transformation with linear regularization are 
graphed. Both spectra corresponding to single species 
(rhodamine 6G or tetramethylrhodamine) have a single 
peak, but the peaks are centered at different values of 
specif ic.brightness.~Jhe peak/of rhodamine 6G is situ- 
ated at about 108 kHz/molecule, whereas the peak of 
tetramethylrhodamine is centered at about 37 kHz/mol- 
ecule. This indicates that a rhodamine 6G molecule is 35 
about 3 times brighter than a tetramethylrhodamine 
molecule. The spectrum corresponding to the mixture of 
the two species has two peaks centered indeed near 
the values obtained for the two species separately. 

Fig. 7 illustrates the residuals corresponding to the ' 40 
measurements of the mixture of rhodamine 6G and 
tetramethylrhodamine. Different methods of data 
processing yield slightly different fit curves (and different 
residuals). The upper graph corresponds to the spec- 
trum of specific brightness obtained by inverse transfor- 45 
mation with linear regularization. The middle graph 
corresponds to the fit curve obtained assuming two spe- 
cies. These two graphs are nearly identical. The experi- 
mentally determined distribution of the number of 
photon counts can formally be fitted under the wrong so 
assumption of single species, which is shown in the 
lower graph, but the fit curve is obviously apart from the 
experimental one. 

Claims ,55 

1. A method for- analyzing samples by measuring 
numbers of photon counts per defined time interval 



in a repetitive mode from light emitted, scattered 
and/or reflected by particles in said sample, and 
determining the distribution of the number of pho- 
ton counts per said time intervals, characterized in 
that the distribution of specific brightnesses of said 
particles is determined from said distribution of the 
number of photon counts. 

. A method according to claim 1 , wherein said parti- 
cles are luminescently labelled or unlabelled mole- 
cules or macromolecules, or dye molecules, 
molecular aggregates, complexes, vesicles, cells, 
viruses, bacteria, beads, or mixtures thereof. 

15 3. A method according to claim 1 and/or 2, wherein 
said particles each carry a number of binding sites 
for luminescent particles. 

4. A method according to at least one of the claims 1 
20 to 3, wherein the measurement volume is only a 

part of the total volume of the sample and said par- 
ticles are diffusing and/or being actively transported 
into and out of said measurement volume and/or 
the sample is actively transported and/or optically 
25 scanned. 

5. A method according to at least one of the claims 1 
to 4, wherein the parameters of the spatial bright- 
ness function characteristic for the optical set-up 
are determined by measuring numbers of photon 
counts per defined time intervals in a repetitive 
mode from light emitted, scattered and/or reflected 

: ~ by a single type of particles ( 

6. A method according to at least one of the claims 1 
to 5, wherein at least one high-aperture microscope 
objective is used in a confocal manner for both 
focussing the incident laser beam and collecting 
light emitted, scattered and/or reflected by particles 
in said sample. 

7. A method according to at least one of the claims 1 
to 6, wherein the size of the pinhole positioned in 
the focal plane of the microscope is used as a mod- 
elling parameter of the spatial brightness function. 

8. A method according to at least one of the claims 1 
to 7, wherein the convergence angle of the incident 
laser beam is used as a modelling parameter of the 
spatial brightness function. 

9. A method according to at least one of the claims 1 
to 8, wherein multiple photon excitation is used to 
excite said particles in said sample. 

10. A method according to at least one of the claims 1 
to 9, wherein the measurement volume is restricted 
by the use of near field optical microscopy, or com- 
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binations thereof with conventional microscopy 
optics. 

>■ F 

1 1 . A method according to at least one of the claims 1 

to 10, wherein the concentration and/or specific 5 
brightness of at least one type of said particles is 
determined/ - 1 

j* . " 

12. A method according to at least one of the claims 1 

to 11, wherein said distribution of the number of W 
photon counts is fitted using a priori information on 
the sample. ■ . ; . 

13. A method according to at least one of the claims 1 

to 12, wherein said distribution of the number of . is 
photon counts is processed by applying an inverse 
transformation with linear regularization. 

14. A method according to at least one of the. claims 1 

to .13, wherein the experimental parameters of the 20 
detector, in particular, the dead time and the after- 
pulsing probability of the detector/ are determined 
by measuring numbers of photon counts per 

defined time intervals in a repetitive mode while the 
detector is exposed to light of constant intensity or 25 
high frequency laser pulses. 

15. A method according to at least one of the claims 1 
to 14, wherein the background count rate of the 
equipment is determined, 30 

16. A method according to at least one of the claims; 1 
to 15, wherein the width of said time interval is'ln 
average smaller than the characteristic correlation 

time of light intensity fluctuations. 35 . 

17. A method according to at least one of the claims 1 
to 16, wherein the width of said time interval is 
selected to yield in average more than one, prefera- ■ 
bly 1 to 10, photon counts per said particle. 40 

18. A method according to at least one of the claims 1 
to 17, wherein the concentration of the sample or 
the size of the measuring volume is selected to 
have in average no more than a few, preferably 45 
about one particle per measuring volume. - 

r 1 

19. A. method according to at least one of the claims 1 

to 18, wherein at least one individual particle is sta- 
tistically analyzed in terms of its specific brightness 50 

which may fluctuate or otherwise change. 
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